Kondo Effect and Josephson Current through a Quantum Dot between Two 

Superconductors 
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We investigate the supercurrent through a quantum dot for the whole range of couplings using 
the numerical renormalization group method. We find that the Josephson current switches abruptly 
from a n- to a 0-phase as the coupling increases. At intermediate couplings the total spin in the 
ground state depends on the phase difference between the two superconductors. Our numerical 
results can explain the crossover in the conductance observed experimentally by Buitelaar et al. 
[Phys. Rev. Lett. 89, 256 801 (2002)]. 

PACS numbers: 74.50.+r, 72.15.Qm, 75.20.Hr 



Introduction. The Kondo effect and superconductiv- 
ity are two of the most extensively studied phenomena 
in condensed matter physics ever since the pioneering 
works by Kondo Q and by Bardeen, Cooper and Schri- 
effer 0] , respectively. When a localized spin is coupled to 
superconducting electrons, the two effects are intermin- 
gled and even richer physics will emerge. The physically 
interesting questions are: Would the Kondo effect sur- 
vive, overcoming the spin-singlet pairing of electrons in 
superconductors (SCs) and the superconducting gap at 
the Fermi level? If it does, how would such a strongly cor- 
related state affect the transport, especially the Joseph- 
son current, between two superconductors? 

The Josephson effect through a strongly interacting re- 
gion with a localized spin was discussed long before by 
Shiba and Soda Q and Glazman and Matveev 0] and 
further elucidated by Spivak and Kivelson 0. The large 
on-site interaction only allows the electrons in a Cooper 
pair to tunnel one by one via virtual processes in which 
the spin ordering of the Cooper pair is reversed, lead- 
ing to a negative Josephson coupling (i.e., a 7r-junction). 
This argument, however, is based on a perturbative idea 
and holds true only for sufficiently weak tunneling. It was 
suggested |^| that as the tunneling increases, the Kondo 
effect produces a collective resonance at the Fermi level. 
As a result, the Josephson current is enhanced by the 
Coulomb repulsion. Moreover, the Josephson coupling is 
expected to be positive (i.e., a 0-junction) since the local- 
ized spin is screened due to the Kondo effect. Based on 
this, Glazman and Matveev Q assumed a strong coupling 
fixed point and derived the Josephson current as a func- 
tion of phase difference. Recently, several approxima- 
tion methods have been used to investigate the transition 
from the 0- to 7r-iunction as a function of the tunneling 
strength Id, LD, IM 1211 : A modified Hartree-Fock approxima- 



Kondo effect. 

In this work, we use a numerical renormalization group 
(NRG) method to investigate thoroughly the 0-n tran- 
sition as well as to examine the argument above sug- 
gested by Glazman and Matveev p| . Based on the NRG 
method, we calculate quantitatively the local properties 
(i.e., the pairing correlation and the single-particle ex- 
citation spectrum) of the quantum dot (QD), the total 
spin in the ground-state wave function, and the Joseph- 
son current as a function of phase difference. Finally, 
we show that our numerical results can explain the ex- 
perimentally observed crossover of the conductance in 
SC-carbon nanotube-SC junctions 

Model. The system consists of a QD with an odd 
number of electrons coupled to two superconducting 
leads (L and R). The study of Kondo effect in such a 
mesoscopic system has recently attracted much interest 
due to its tunability. As alr eady demonstrated exper- 
imentally with normal leads [ll|, it allows for various 
tests of Kondo physics, which are difficult in bulk solids. 
The two leads are regarded to be standard s-wave super- 
conductors (SCs) and described by the BCS Hamiltonian 



^bcs = ^2 ^2 Q,kc[k,a C «,k CT 
£=L,R k,<r 

-EE( A ^ c U, T 4-ka + H ' (1) 

l k 

where c\ k a (c^k i<T ) creates (destroys) an electron with 
energy e^k, momentum ?ik, and spin a on the lead I. Ai 
is the superconducting gap and fa is the phase of the 
superconducting order parameter. The QD is described 
by an Anderson-type impurity model 



(2) 



tion |fjj, a non-crossing approximation (NCA) jZJ, and a 
variational method || predict a 0-7T transition, whereas 
the slave-boson mean-field theory || always favors the 



which is widely adopted for sufficiently small quantum 
dots. In Eq. J5J d\ and d a are electron creation and 



annihilation operators on the QD. The level position e^, 
measured from the Fermi energy Ep of the two leads 
(throughout the paper every energy is measured from 
Ep), can be tuned by an external gate voltage. The 
interaction U is order of charging energy e 2 /2C (C is the 
capacitance of the QD). The coupling between the QD 
and the SCs is described by the tunneling Hamiltonian 



(3) 



k.a 



Putting all together the Hamiltonian for the whole sys- 
tem is given by 3f — J#qd + J^hcs + ■ 

We take a few simplifications to make clearer the 
physical interpretation of the results below. The two 
SCs are assumed to be identical (fL.k = e_R,k = £k 
and Al — An = A) except for a finite phase differ- 
ence <j) = 4>l — 4>r\ without loss of generality we put 
<f>L = —<f>R = <A/2. In the normal state, the conduction 
bands on the leads are symmetric with a flat density of 
states No and the width D above and below the Fermi 
energy. We also put e<z — —U/2 in J^qn, Eq. (J2J; it has 
been checked that an asymmetric model (ed ^ —U/2) 
gives the qualitatively same results for physical quan- 
tities of our concern. We only consider the symmetric 
junction, Vl = Vr = V. The coupling to the leads is 
well characterized by the single parameter T = 2ttNoV 2 . 
Below we will distinguish the strong (Tk A) and the 
weak (Tk <C A) coupling limits by the ratio between 
the superconducting gap A and the normal-state Kondo 
temperature Tk (fcs = 1) given by ^3] 



7T 1 H 

2r V u 



(4) 



Following the standard NRG procedures 0, ex- 
tended to superconducting leads [l^, we evaluate the 
various physical quantities from the recursion relation 

N a jp,,N+l,(T + h.C. J 

A w / 2 A, (fl N+u fl, N+ld + h.c) (5) 
with the initial Hamiltonian given by 



Mo = —j= 
VA 



fj,=e,o a 



(6) 



Here the fermion operators fij.,N,a have been introduced 
as a result of the logarithmic discretization and the 
accompanying canonical transformation, A is the log- 
arithmic discretization parameter (we choose A = 2), 
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FIG. 1: The pairing correlation on the quantum dot, Ad = 
(djdj), as a function of A/Tk- Inset: plot of bare (not nor- 
malized) values of Ad. We have chosen ed = —U/2 = —0.1D 
and T = 0.04D. 



&v ~ 1 0, and 



V R 




(7) 



sin(<£/4) . 



with ( = 1+ 2 / A . The Hamiltonians Jtfjy in Eq. (JSJ 
have been rescaled for numerical accuracy. The original 
Hamiltonian is recovered by Jif / D = limjv— ><x> -J^n/^n 

with y N = ca^- 1 )/ 2 . 

Proximity effect. To see how superconductivity on 
the leads affects the interacting QD in the strong and 
weak coupling limits, we first examine the local proper- 
ties on the QD with zero phase difference (<f> = 0) |16| . 
Figure shows the local pair correlation A d = (dtd\) as 
a function of A/Tk- As expected, the local pair correla- 
tion Ad vanishes with A, and gets smaller (even vanishes 
when U — > oo) as A — > oo; see Fig.^(inset). An interest- 
ing aspect of Ad is the sign change at A = A c ~ 2.4T^, 
which suggests that the physical properties are different 
in the strong (Tk ^> A) and the weak (Tk *C A) cou- 
pling limits. Indeed we see (from the NRG calculation) 
that the ground-state wave function of the whole sys- 
tem is of spin singlet (the localized spin is screened out) 
for A < A c and of spin doublet (the SCs form Cooper 
pairs separately and the localized spin is left unscreened) 
for A > A c . The negative sign in Ad in the weak cou- 
pling limit can be explained by a simple second-order 
perturbation theory, while the positive on in the strong- 
coupling limit is expected when there is a resonance chan- 
nel for Cooper-pair tunneling Therefore, it seems 
quite plausible to argue that in the strong coupling limit 
the Kondo resonance develops even in the presence of 
the superconducting gap in the conduction band and the 
proximity effect arises through the resonance; see also 
the discussion of the Josephson current below. Putting 
it another way, the local moment of spin 1/2 induces a 
negative Ad for weak couplings, but as the coupling in- 
creases it is screened and a positive A^ is recovered. 
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FIG. 2: The single-particle excitation spectrum on the quan- 
tum dot. e d = -U/2 = -0.LD and T = 0.04D (T K = 
0.0089D). 



This interpretation is further supported by the single- 
particle excitation spectra Ad(E) on the QD, as shown in 
Fig.Elfor different values of A/T K . In Fig. 0(a) A(E) for 
zero phase difference 0] is shown and we observe a qual- 
itative change of the spectrum when A becomes smaller 
than Tk- A localized state below the superconducting 
gap appears for A > Tk, whereas the spectrum has a gap 
of the order of A in the other limit. The other panels in 
Fig. [21 show the phase-dependent density of states in the 
sub-gap regime. We clearly observe a phase-dependent 
formation of an Andreev bound state. For A/Tk = 0.1 
the Andreev state emerges from the gap with increasing 
phase and reaches the smallest energy for <f> = 7T, which is 
reminiscent of a usual superconducting junction. In weak 
coupling limit, A = IOTk, we observe an opposite phase- 
dependence, which is similar to the predicted 7r-junction 
behavior p|. For an intermediate coupling, A/Tr- = 1.8, 
there is always a localized state below the gap, which has 
a non-monotonic phase-dependence. In the following, we 
will discuss the Josephson current through the quantum 
dot. 

Josephson current. We now turn to the Josephson 
current through the QD in the presence of a finite phase 
difference 4>. Within the NRG method, the Josephson 
current can be conveniently calculated by the relation [lij 



J short 



DT 

2ttA 2 



[sin(0/4)J e +cos(</>/4)J o ] (8) 



short — 



with J M = J2 a ( d ifvfl,a + h.c.) (jjl = e,o). Here J® 
eA/h is the critical current of a transparent single- mode 
junction 

Figure |3 shows the Josephson current as a function 
of phase difference 4> between the two superconducting 
leads for different values of ratio A/Tr-. In the weak 
coupling limit (Tk <C A), it is clearly seen from Fig. [3] (a) 
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FIG. 3: Josephson current I s (4>) (in units of J= hort = eA/h) 
as a function of phase different <j> (a) for A/Tk = 10 and (b) 
for A/T K = 0.1. (c) Same curves for A /Tk = 1-6, 1.8, 2.0, 
and 2.2 (near the 0-7r junction transition point), (d) Critical 
current in the Kondo regime. We put ed = —U/2 — —0.1D 
and r = 0.04D. Inset: conductance resulting from the RSJ- 
model (see text). 



that the effective Josephson coupling is negative (i.e., a 
7r-junction) 0, 0, 0J LI ■ In addition, the supercurrent- 
phase relation is very close to a sinusoidal function, like 
typical "tunneling junctions" |17| . We also report that 
the ground state is a doublet for any phase difference <fi. 

In the strong coupling limit (Tk S> A), on the other 
hand, the Josephson coupling is positive 0, 0; see 
Fig- El (b)- Another remarkable thing is that the current- 
phase relation is highly non-sinusoidal and reminiscent of 
the current-phase relation in the short junction limit |l7| . 
Furthermore, the critical current approaches the unitary 
limit i^ hort of "short junctions" ^3 as the coupling grows 
stronger (A/Tk — ► 0), as shown in Fig. 01(d). These re- 
sults suggest again that in the strong coupling limit the 
Kondo resonance develops at the Fermi level and Cooper 
pairs tunnel resonantly through it. Naturally, the ground 
state turns out to be a spin singlet for any </>. It should 
be stressed here that although the Kondo effect manifests 
itself as a resonance channel for the Cooper-pair tunnel- 
ing, the Kondo peak of width Tk in the quasi-particlc 
excitation spectrum is suppressed (showing a gap) below 
the energy scale of order A (<§; Tk)] see Fig- EI 

Another interesting regime is the intermediate one 
(A - T K )- As demonstrated in Fig.El(c), for A - T K the 
curve of Xs(</>) breaks into three distinct segments. The 
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central segment resembles that of a ballistic short junc- 
tion ^3 1 while the two surrounding segments are parts of 
a 7r-junction curve 0. Namely, the critical value AJ6) 
depends on <j> with A c (<f>) > A c (<f>') for \d)\ < |</>'| \2£; 
for example, A c (0.3tt) w 1.6 and A c (0) re 2.4. Evidently, 
the NRG results show that the ground state is a spin sin- 
glet in the central segments (A < A c ((f))) and a doublet 
in the other (A > A c ). 

Experiments. In the experiments of Buitelaar et al. 
|l0| the interplay between superconductivity and Kondo 
physics was observed in non-e quili brium transport (mul- 
tiple Andreev reflections) HSHil, but no supercurrent 
was measured. However, the absence of a dissipationlcss 
branch in the IV is not surprising in such (intrinsically) 
small junctions. Indeed thermal or quantum fluctuations 
in connection with a resistive environment can lead to 
a finite resistance [24j. In Ref. the "quality factor" 
i?ArC(2e/ c /?iC) 1 / 2 , governing the dynamics of the cor- 
responding resistively-shunted junction (RSJ) model, is 
always smaller than 1 [2j| and the junction is therefore 
overdamped. In this limit the measured resistance Gs 
is directly related to the current-phase relation, roughly 
like Gs/Gn ~ cxp(TiI c /eT) 26] . This enables us to re- 
late our results of Fig. to the measured crossover of 
the conductance as function of Tk/A, see Fig. 4 of flflj . 
For the experimental temperature T — 50mK and gap 
parameter A ~ 1.2K, the calculated critical current in 
Fig. El (d) means that the factor hI c /eT becomes much 
larger than 1 in the Kondo regime A <C Tk, when the 
transparent junction limit is reached. Thus, the experi- 
mentally observed crossover to Gs > Gn in this limit is 
a manifestation of the supercurrent approaching the uni- 
tary limit eA/h. The inset of Fig. [3] (d) shows the con- 
ductance as a function of Tk/A and that the crossover 
appears for Tk ~ 0.5A, which is in quite good agreement 
with the experimental result of Ref. ^(J ■ 

Conclusion. We have studied the Josephson current 
and the proximity effect on the QD coupled to two SCs 
in a whole range of coupling. Our results exhibit a tran- 
sition from the weak to the strong coupling limit, which 
occurs when A ~ Tk- In the weak coupling limit, su- 
perconductivity dominates the Kondo physics, and the 
tunneling of Cooper pairs can be treated perturbatively. 
The system is a 7r-junction, the pairing correlation on the 
QD is negative, and the ground state is a spin doublet. In 
the strong coupling limit, the Kondo effect becomes im- 
portant and manifests itself as a resonance channel for the 
Cooper-pair tunneling. This leads to a positive Joseph- 
son coupling (0-junction) and positive pairing correlation 
on the QD. Here the Kondo effect in the presence of su- 
perconductivity is distinguished from the usual one with 
normal leads in that the Kondo peak in the quasi-particle 
excitation spectrum is suppressed completely (exhibiting 
a gap) for energies below the superconducting gap. 
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